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摘 要 为 了 探究 不 同 转速 对 某 亚 音 轴 流 压气 机 叶 尖 泄漏 流 非 定常 性 的 影响 , 采用 动态 压力 测量 技术 获得 不 同 转速 机 匣 
壁 处 压力 信号 . 试验 测量 表明 , 在 近 失速 工 况 压力 信号 的 频谱 出 现 表 征 旋转 不 稳定 性 的 特征 频率 驼峰 . 其 驼峰 峰值 处 频率 
随 转速 的 提高 相对 减 小 . 在 此 基础 ， 分 别 进行 不 同 转速 的 多 通道 数值 模拟 . 数值 模拟 中 ,， 叶 尖 区 的 静 压 监测 信号 频谱 与 试 
验 具有 一 致 性 , 转子 叶 尖 存在 逆转 子 传播 的 周 向 行 波 . 该 周 向 行 波 对 应 了 试验 中 的 旋转 不 稳定 现象 . 进一步 的 数值 模拟 流 
场 分 析 表 明 , 该 周 向 行 波 在 53 史 设计 转速 时 是 源 于 叶 顶 泄漏 流 自我 维持 的 非 定常 性 . 在 71% 设 计 转 速 时 , 在 某 些 时 刻 泄漏 
涡 发 生 破 碎 并 参与 到 周 向 行 波 的 形成 过 程 . 
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. Abstract In order to explore yeaee of the different rotor speed on the unsteady flow near tip 
es region in a subsonic axial 了 ipressor rotor, the casing wall static pressure signals were measured 
ER by the dynamical presswae tkansducers at 53% and 71% design rotor speeds respectively. In the 
pressure frequency OA were to be found as characteristic humps below the blade passing fre- 
quency near stall at different rotor speeds. The peak frequency of the Characteristic hump shifts lower 


with the increase of the rotor speed. Based on the experimental Basurement, the multi-passage sim- 
ulations were conducted at different rotor speeds. The charactéyistic frequencies of the static pressure 


oscillations near tip region at different rotor speeds were gaitent with that in measurements. It can 
be found that circumferential traveling waves rotat Agaifst the rotor near tip region at the two rotor 
speeds. The circumferential traveling waves can be rd as RI. The details of the flow fields indicated 
that RI originated from the self-induced of 半 i eakage flow at 53% design rotor speed. However, 
the formation of RI was accompanied by Yeakdown of tip leakage vortex at 71% design rotor speed. 
Key words subsonic compressor; axial flow; unsteady flow; tip leakage flow; experimental measure- 


ment; numerical simulation:; 


0 引 言 


研究 表明 ， 压 气 机 叶 尖 区 的 非 定常 流动 是 导致 
压气 机 失 稳 的 主要 因素 之 一 。 科研 人 员 对 此 进行 了 
大 量 的 试验 及 数值 模拟 研究 。 早期 针对 失速 起 始 工 
况 的 研究 ， 学 者 们 分 别提 出 了 两 种 失速 先兆 波 的 概 
念 ， 既 模 态 波 站 和 突 尖 波 站 ,对 应 突 尖 波 型 失速 
先兆 的 机 理 研究 ， 后续 的 学 者 发 现 这 与 泄漏 流 或 湛 
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漏 涡 的 运动 密切 相关 中 .这 期 间 Voig 就 提出 了 判 
断 突 尖 波 型 失速 的 准则 ， 

而 近年 来 ， 在 压气 机 靠近 失速 边界 的 稳定 工 
况 ， 对 间 隐 泄漏 流 非 定常 性 的 研究 在 国内 外 也 得 
到 了 关注 。Furukawa 等 [" 引 分 别 采 用 定常 数值 模 
拟 和 非 定常 数值 模拟 对 某 低速 轴 流 压气 机 的 间隙 泄 
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漏 涡 破 雁 现 象 进行 研究 ， 并 用 其 提出 的 潍 漏 涡 破 碎 
理论 解释 了 其 流 场 叶 尖 区 的 非 定 常 现象 以 及 分 析 泄 
漏 涡 破碎 对 压气 机 性 能 的 影响 .Hoffman 等 四 以 
及 Schlechtriemllol 相应 的 研究 将 泄漏 涡 破碎 与 压 
气 机 失 稳 之 间 联 系 起 来 。 此 外 ，Zhang 等 1 针 
对 高 载荷 的 低速 转子 的 研究 认为 间 际 泄漏 流 非 定 
常 性 是 主 来 流 与 泄漏 流 相互 作 用 所 引起 的 称 为 自 诱 
导 型 非 定常 性 导致 的 ， 并 得 到 Tong 等 03 在 试 
验 上 的 证 明 。 Duts.14 等 将 泄漏 流 自 诱 导 非 定常 性 
进一步 拓展 到 跨 音 速 条 件 下 ， 其 单 通道 非 定常 数 
值 模拟 的 近 失 速 工 况 点 均 发 现 了 泄漏 流 自 诱导 非 定 
常 性 ， 

在 对 近 失 速 稳定 工 况 叶 尖 区 非 定常 性 的 研究 中 ， 
科研 人 员 发 现 了 一 种 类 似 突 尖 波 的 短 尺 度 周 向 扰动 
波 ， 并 称 之 为 旋转 不 稳定 性 (Rotating Instability, 
RD). 在 揭示 RI 产生 机 理 的 研究 中 , Mailach 等 [15,19] 
认为 RI 源 于 叶 尖 泄漏 涡 的 波动 . 而 Mirz 等 上 7 针 
某 轴 流通 风机 的 研究 则 认为 RI 源 于 轴 向 反 流 和 来 
流 在 相互 作用 形成 一 个 新 径 向 涡 的 周期 性 运动 . 


常数 值 模拟 和 非 定常 数值 模拟 对 叶 尖 区 全 
深入 的 研究 , 并 得 到 与 Mirz 等 0 


Dh 乎 表明 了 RI 
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与 压气 机 的 突 尖 波 失速 先兆 的 关联 . 

我 们 发 现 , 与 之 前 所 提 Qw 等 049 相同 的 研究 对 
象 中 , Biela 等 Po 的 试验 研究 也 发 现 RI 现象 . 这 说 
明了 Du 等 04 所 观察 到 的 自我 诱导 非 定常 性 与 RI 
可 能 存在 密切 关联 . 但 Du 的 研究 仅 在 单 通道 的 数 
值 模拟 下 进行 , 无 法 解释 该 现象 和 RI 的 关系 . 因此 
本 文 以 某 亚 音 轴 流 压气 机 孤立 转子 为 研究 对 象 ， 进 


行 两 个 转速 (53% 设 计 转 速 和 71% 设 计 转 速 ) 的 试验 < . 


测 ,并 在 此 基础 进行 多 通道 的 数值 模拟 . 其 目的 xc 
是 揭示 自我 诱导 非 定常 性 与 RI 的 关联 , 二 是 扰 完 不 
同 转速 对 RI 产生 过 程 的 影响 ， 


1 研究 对 象 及 试验 测量 


1.1 研究 对 象 

本 文 研究 对 象 为 西北 工业 大 学 亚 音 轴 流 压 气 机 
孤立 转子 试验 台 , 该 试验 台 设 计 转 速 为 15200 r/min， 
轮 载 比 为 0.6, 叶 尖 直径 为 298 mm 以 及 叶 尖 稠度 为 
0.961, 其 详细 的 设计 参数 和 几何 参数 见 文献 [21], 图 
1 给 出 了 该 试验 台 的 结构 简 图 . 本 文 试验 测试 是 在 
孤立 转子 环境 下 进行 的 ， 即 把 静 子 后 移 至 排 气 管道 
进口 的 位 置 . 试验 的 叶 顶 径 向 间 际 为 0.35 mm, 其 转 
速 为 设计 转速 的 53% 和 71%%。 
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1. Motor 2.Accerlator 3.Torque meter 4. Measurement plane 
5. Rotor 6. Stator 


图 1 轴 流 压气 机 试验 台 结构 简 图 


Fig. 1 _ Cross-sectional diagram of the test Tig 


7. Discharge plenum 8.Throttle valve 


试验 时 调节 试验 台 出 口 节气 门 的 开 度 可 以 控制 
通过 试验 转子 的 流量 , 流量 的 测量 由 安装 在 排 气 管 
道中 的 孔 板式 流量 计 完 成 . 气流 通过 转子 的 压 比 由 
安装 在 进 、 出 口 测量 截面 上 的 气动 探 针 获得 . 通过 
关 小 节气 门 提高 压气 机 的 出 口 背 压 ， 从 而 改变 压气 


机 运行 同时 在 相应 工 况 处 进行 了 稳 态 数据 采 
集 给 转 子 的 稳 态 特性 和 稳定 工作 的 流量 范 


洛 稳 态 测量 的 基础 上 ， 从 峰值 效率 工 况 (PE) 微 
前 节气 门 直至 失速 边界 ， 在 每 个 稳定 工 况 采集 动态 
静 压 信号 , 采样 频率 为 100 kHz, 每 个 工 况 的 数据 存 
储 时 间 为 3 s. 

其 中 机 面壁 上 的 动态 压力 测量 是 由 高 频 响 动 
态 压力 测量 系统 完成 。 它 由 高 频 动态 压力 传 感 
器 (Kulite XCQ-080-5G)、 信号 调理 器 (Preston 
8300AU)、 信 号 采集 、 分 析 模 块 (LMS) 以 及 计算 
机 组 成 . 动态 压 xX 传 感 器 重 直 安装 于 正 对 转子 的 机 
匣 壁 孔 内 ， 基 安 半 位 置 示意 见 图 2。 以 叶 顶 叶片 
前 缘 为 基 8V 测 点 1、2、3、4、5、6、7 分 别 为 
一 25 XY%, 19.4%，38.8%，58.2%，77.6% 和 97% 叶 


图 2 高 频 动态 传感器 放置 位 置 图 


Fig. 2 Measured and predicted total pressure ratio 


试验 通过 关 小 节气 门 提高 压气 机 出 口 背 压 改变 
压气 机 的 运行 工 况 ,， 并 在 相应 的 工 况 点 处 进行 了 稳 
态 数据 采集 ， 得 到 试验 转子 稳 态 特性 和 稳定 工作 的 
流量 范围 .在 稳 态 测量 的 基础 上 ， 从 峰值 效率 工 况 
(标记 为 PE) 开始 进行 机 匣 壁 动态 压力 采集 , 通过 微 
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调节 气门 位 置 改 变 工 况 ， 压 气 机 进入 旋转 失速 的 前 
的 最 后 一 个 稳定 工 况 定 为 实验 边界 工 况 ， 该 工 况 点 
在 后 续 的 转子 性 能 图 (图 5) 中 以 NS 标识 .动态 压 
力 采样 频率 为 100 kHz， 每 个 工 况 的 数据 存储 时 间 
为 3 s。 
2 数值 模拟 方法 
2.1 计算 方法 

本 文 数值 模拟 采用 NUMECA Fine/Turbo 模 
块 , 采用 Spalart-Allmaras 方程 来 封闭 时 均 Navior- 
Stokes 方程 。 有限 体积 中 心 格式 对 控制 方程 进行 离 
散 . 定常 计算 采用 显 式 4 阶 Runge-Kutta 方法 . 非 
定常 计算 采用 隐 式 双 时 间 步 方法 ， 其 中 跟 据 计 算 结 
果 无 依赖 性 的 临界 时 间 步 长 为 1.23x10-5, 叶片 通过 
一 个 栅 距 给 定 20 个 物理 时 间 步 , 每 一 物理 时 间 步 下 
的 虚拟 时 间 步 为 50. 计算 空间 离散 格式 采用 基于 通 
量 差分 离散 法 的 二 阶 迎 风格 式 ， 并 采用 Van albada 
限制 器 限制 解 的 振荡 。 另 外 采用 多 重 网 格 法 、 当 地 
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时 间 步 长 法 及 险 式 开关 光大 竺 方法 加 问 收 但 训 克 


算 效率 ， 次 
2.2 计算 网 格 

各 考 记 新 要 站 针对 次 让 娄 开 的 由 
密度 对 的 研究 , 本文 单 通道 网 格 总 前 溃 点 数 超过 61 
万 满足 网 格 无 关 性 要 求 ， 图 C3 各 了 本 文 计算 所 
采用 的 五 通道 网 格 . 与 叶片 扯 交 的 转动 块 采 用 HI 型 
网 格 拓扑 结构 转动 块 的 前 面 为 型 静止 块 , 设置 
萝 止 块 的 目的 是 为 了 实现 静止 坐标 系 下 压力 监测 点 
的 布置 . 为 了 便于 后 续 分 析 , 如 图 3 将 叶片 编号 , 其 
中 叶片 1 和 叶片 2 之 间 的 通道 称 为 通道 L2, 并 以 此 
类 推 


R/S Interface 


图 3 计算 网 格 拓扑 结构 
Fig. 3 Computational grid 


2.3 数值 探 针 布置 
非 定常 数值 计算 时 ， 如 图 4 给 出 了 在 叶 尖 区 
即 98.15% 通 道 高 的 截面 轴 向 数值 探 针 的 布置 图 . 为 


等 : 亚 音 轴 流 压气 机 转子 转速 对 叶 尖 区 非 定常 流 的 影响 


了 与 试验 中 动态 压力 测量 点 位 置 相 对 应 ， 在 转动 
块 通道 3-4 内 近 压 力 面 和 近 吸 力 面 沿 流向 分 别 在 
0%，19.4%%，38.8%，58.2%，77.6% 和 97% 叶 顶 轴 向 
弦 长 处 布置 数值 探 针 (对 应 试验 测 点 2-7), 在 静止 块 
正 对 通道 3-4 的 -25% 叶 顶 轴 向 弦 长 处 (对 应 试验 测 
点 1) 以 及 其 他 通道 位 置 38.8% 轴 向 弦 长 处 也 进行 数 
值 探 针 布置 ， 


、Bl B2 B3 B4 B5 
i Em ge en i 28.8%Cx 


—25%C, R/S Interface 
图 4 数值 探 针 布置 示意 图 


Fig. 4 The arrangement of numerical probes on the cross 


surface at blade tip 


2.4 读 
六 值 模拟 分 别 在 两 个 转速 下 进行 ,其 计算 步 又 
阁 行 定常 计算 获得 压气 机 稳定 工作 特性 线 , 从 定 

计算 的 峰值 效率 工 况 点 开始 进行 非 定 常 计算 ， 通 
过 微调 计算 域 出 口 背 压 的 方式 台 近 非 定常 计算 数值 
边界 点 。 
3 试验 及 数值 模拟 结果 分 析 
3.1 总 性 能 对 比 

图 5 比较 了 两 个 转速 的 计算 和 试验 压 比 性 能 图 . 
图 中 的 非 定常 计算 绍 果 为 一 个 计算 周期 内 的 时 均值 . 
从 图 中 可 以 明显 地 看 出 , 在 53% 设 计 转 速 和 71% 设 
计 转 速 的 数 稚 模拟 结果 与 试验 数据 变化 趋势 都 能 较 
好 也 “其 中 PE 为 峰值 效率 工 况 ， NS 为 失速 边 

饮 . 在 对 试验 及 数值 模拟 各 个 工 况 测 得 动态 压 
信和 号 进行 频谱 分 析 的 基础 上 ， 可 以 得 到 如 下 两 个 
结论 : 一 是 在 两 个 不 同 转速 下 都 是 在 靠近 失速 边界 
极 小 的 流量 范围 的 工 况 出 现 非 定常 性 ; 二 是 随 着 转 
速 的 提高 出 现 非 定 常 性 工 况 的 流量 范围 变 大 . 
3.2 试验 测量 信号 分 析 

图 6 给 出 了 两 个 转速 下 各 自 PE 工 况 和 NS 工 
况 在 探 针 4 测 得 的 动态 压力 信号 功率 谱 分 析 比 较 图 . 
图 中 的 主要 的 特征 频率 是 叶片 通过 频率 (Blade Pass- 
ing Frequency, BPF) 及 其 倍 频 . 而 相 比 PE 工 况 , 在 
NS 工 况 无 论 是 53% 设 计 转 速 还 是 71% 设 计 转 速 , 探 
针 4 测 得 的 动态 压力 信号 频谱 上 都 出 现 了 一 个 低频 
的 驼峰 状 的 特征 频率 ， 其 驼峰 峰值 处 频率 分 别 约 为 
0.47 倍 BPF 和 0.34BPF。 综 上 可 知 流 场 中 出 现 的 
非 定 常 性 了 是 RIn7, 同时 也 说 明 随 着 转速 的 提高 ， 
在 失速 边界 出 现 的 非 定 常 性 其 特征 频率 相对 性 的 
减 小 。 
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Total pressure ratio 
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(a) 53% 设 计 转 速 (b) 71% design rotor speed 
(a) 53% design rotor speed 图 6 探 针 4 测量 信号 功率 谱 分 析 


1.13 rey 和 了 Fig. ectra of static pressure signals detecte e 
: : ] ' l a ee 
“各 数值 模拟 信号 分 析 
on 两 个 转速 的 数值 模拟 流 场 中 出 现 非 定 常 性 的 
A 流量 范围 分 别 为 (Mass flow==2.91~3.00 kg.s-1) 和 
(Mass flow 二 2.94~4.16 kg:s-!), 见 图 5。 为 了 与 试验 
作 比 较 , 分 别 选 取 两 个 转速 数值 模拟 的 失速 边界 点 
(Mass flow=3.00 kg:s-!) 和 (Mass flow=4.16 kg:s-!) 
作为 数值 分 析 对 象 . 图 7 给 出 了 两 个 转速 相对 坐标 


Total pressure ratio 
一 
i 
© 
(LA 


1.085E 系 下 38.8% 轴 向 弦 长 处 数值 探 针 的 信号 FFT 频谱 分 
1 析 图 ,从 图 7 可 允许 出 ， 两 个 转速 的 数值 模拟 流 志 
Mass flow/kg:s”! 中 存在 非 定 JE， 结合 图 8 中 静 压 监测 点 的 压力 波 

(b) 71% 设 计 转 束 动 周期 (前 者 了 = 61 后 者 了 = 78), 说 明 相 对 坐标 系 
Ta 71% 设 计 转 速 时 比 53% 设 计 转 束 

图 5 实验 和 计算 的 总 压 比 特性 泣 长 ， 既 静 压 扰动 频率 随 着 转速 提高 而 减 小 ,这 和 

人 试验 存在 着 一 致 性 . 此 外 由 图 8 也 可 以 看 到 , 在 两 种 


Fig. 5 Measured and predicted total pressure ratio SS 


” 转速 中 相 邻 通道 的 静 压 波动 存在 一 个 固定 的 时 间 相 
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图 7 相对 坐标 系 98.15% 通 道 高 度 近 压力 面 侧 38.8% 轴 向 弦 长 


Fig. 7 Spe 


106500 


处 数值 监测 信号 频谱 分 析 


ctra of static pressure signals detected at 98.15% 
span and 38.8% tip axial chord 
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(a) 53% 设 计 转 速 

(a) 53% design rotor speed 
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(b) 71% 设 计 转 速 

(b) 71% design rotor speed 


图 8 相对 多 


E 标 系 98.15% 通 道 高 近 压 力 面 侧 38.8%% 轴 向 弦 长 处 


5 个 通道 测 点 的 静 压 监测 结果 


Fig. 8 Pressure oscillations taken 位 om probes in 5 passages 


at 98.15% span and 38.8% tip axial chord near pressure side 


位 差 . 这 说 明 两 种 转速 近 失 速 工 况 的 流 场 中 ， 转 子 
叶 尖 区 域 均 出 现 了 周 向 传播 的 扰动 行 波 ， 

试验 是 在 机 匣 壁 上 进行 的 动态 压力 信号 测量 ， 
压力 探 针 处 于 静止 坐标 系 下 . 为 了 与 试验 进一步 比 
较 , 图 9 给 出 数值 模拟 在 绝对 坐标 系 下 的 监测 点 信 
号 FFT 频谱 分 析 图 . 两 个 转速 数值 模拟 绝对 坐标 系 
下 , 流 场 中 非 定常 性 的 频率 大 小 分 别 为 0.47 倍 BPF 
和 0.34 倍 BPF, 这 与 试验 的 信号 分 析 频 谱 图 中 驼峰 
ee 综 上 分 析 ， 即 试验 测量 信 
号 、 数 值 监测 信号 频谱 分 析 以 及 数值 模拟 中 叶 尖 区 
nn 
的 数值 模拟 成 功 的 模拟 出 了 试验 中 的 RI 非 定常 性 ， 
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图 9 绝对 坐标 系数 值 监测 信号 频谱 分 析 


Fig. 9 Spectra of static pressure signals detected at 


stationary part 


3.4 数值 模拟 流 场 分 析 
上 述 分 析 已 经 说 明 数 值 模拟 和 试验 的 一 致 性 ， 
那么 以 下 将 分 别 对 两 个 转速 的 内 部 流 场 进行 详细 的 。 
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图 10 和 图 11 分 别 给 出 了 53% 和 71% 设 计 转 速 数值 
模拟 失速 边界 工 况 五 通道 一 个 扰动 周期 泄漏 流 及 各 
叶片 压力 面 静 压 分 布 图 , 扰动 周期 7 取 为 压力 的 近 
波 谷 到 近 波 谷 , 如 图 8 所 示 . 更 压 系数 表达 式 如 下 : 


OC» = (p — po) /(0.5pU7) (1) 


其 中 2 为 静 压 ; po 为 参考 压力 标准 大 气压 ; p 为 密 
度 ; Ui 为 叶 尖 线 速度 . 
3.4.1 53% 设 计 转 速 近 失 速 工 况 流 场 分 析 

分 析 可 知 主 泄漏 在 近 失 速 边 界 时 与 非 定常 性 未 
出 现 的 工 况 相 比 并 没有 发 生 明 显 的 变化 ， 而 尾 缘 的 
泄漏 流 却 发 生 周 期 性 运动 的 变化 . 因此 为 了 便于 分 
析 将 叶 人 尖 间 际 泄漏 流 分 为 a 和 b 两 段 . a 段 为 0% 到 
40% 轴 向 弦 长 段 ,b 段 为 40% 到 100% 轴 向 弦 长 段 以 段 . 

在 上 = (0/6)T 时 刻 , 通道 3-4 中 叶片 4a 段 的 泄 
漏 流 冲击 到 叶片 3 压力 面 约 38.8 多 轴 向 弦 长 处 后 贴 
着 压力 面 流 出 通道 (其 轨迹 如 图 中 黑色 实 线 所 示 ) 并 
在 约 38.8% 轴 向 豆 长 处 形成 一 个 低压 区 (图 中 由 圆圈 
I 示 出 )。 


两 部 分 ,一 部 分 仍 是 冲击 在 叶片 3 压力 

力 面 流 出 通道 , 另 一 部 分 则 是 发 生 了 汽 漏 . 这 
是 由 于 叶片 4 压力 低压 区 后 的 敲 奈 局 ~( 由 圆圈 工 示 
出 ) 进一步 向 叶片 尾 缘 移动 ， 片 4b 段 处 压 差 
增 大 ,该 处 汇 漏 流 获 得 较 E 量 ， 最 终 使 得 部 分 


a 
在 t= (1/6)T 时 刻 叶片 4b 段 的 泄漏 流 
压 


叶片 4b 段 泄漏 流 在 叶片 3 发 生 第 二 次 泄漏 , 冲击 叶 
片 3 压力 面 的 叶片 4b 段 泄漏 流 强度 随 之 减 小 , 其 产 
生 的 低压 区 也 相应 的 减弱 的 . 

到 t= (2/6)7 时 叶片 4 压力 面 b 段 处 的 压力 进 
一 步 增 大 , 使 得 全 段 的 叶片 4p 段 泄漏 流 在 叶片 3 发 
生 第 二 次 泄漏 ， 叶 片 3 压力 面 低压 区 几乎 消失 . 叶 
片 4b 段 汽 漏 流 发 生 第 二 次 泄漏 后 , 冲击 叶片 2 压力 
面前 缘 并 形成 叶片 2 压力 面 新 的 低压 区 。 随 着 时 间 
的 发 展 ， 由 于 叶片 4 压力 面 低压 区 移动 到 尾 缘 的 缘 
故 , 叶片 4b 段 处 的 压力 开始 减 小 , 从 而 使 得 叶片 4b 
段 泄漏 流 的 强度 逐渐 无 法 使 得 自己 在 叶片 3 发 生 全 
段 的 第 二 次 泄漏 , 在 上 = (3/6)7 时 刻 叶 片 4b 段 泄 漏 
流 开始 变 得 澳 散 , 到 上 = (4/6)7 时 叶片 4b 段 泄漏 
流 “ 退 化 ” 变 回 两 部 分 ， 此 时 叶片 5b 段 泄漏 流 发 生 


其 冲击 叶片 3 压力 面 形成 的 低压 区 
与 叶片 5b 段 泄 漏 流 第 二 次 泄漏 冲击 叶片 3 压力 面 
形成 的 新 低压 区 重合 . 而 随 着 叶片 4 压力 面 新 的 高 
压 区 向 叶片 尾 缘 的 移动 , 叶片 4b 段 泄 漏 流 强度 重新 
增强 . 此 时 原先 t = (0/6)T 时 刻 通 道 3-4 的 泄漏 流 
状态 已 经 传播 到 了 通道 2-3 中 (黑色 虚线 所 示 )。 到 
t = (6/6)7 时 刻 流 场 与 上 = (0/6)7 时 刻 一 致 ， 此 时 
叶片 5b RS 叶片 3 压力 面 

CG 


0.4 (a) 大 (0/6) 了 


(b) (1/0)T 


(FQOT (qd) £3/OT 


(FAOT (ESOT 


图 10 53% 设计 转速 数值 模拟 失速 边界 工 况 一 个 静 压 扰动 周期 近 压力 面 侧 静 压 分 布 及 泄漏 流 线 图 


Fig. 10 Chord-wise distribution of static pressure coefficient on the pressure side at 98.15% span and TLF streamlines in passage 


3-4 during one period of pressure oscillation for NS at 53% design rotor speed 
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图 11 71% 设计 转速 数值 模拟 失速 边 


Fig. 11 Chord-wise distribution of static pressure 


(d) 六 3/5)7 (EAT A/5)T 
个 静 压 扰动 周期 近 压 力 面 侧 静 压 分 布 及 泄漏 流 线 图 


lent on the pressure side at 98.15% span and TLF streamlines in passage 


3-4 during one perio sure oscillation for NS at 71% design rotor speed 


必 


人 
tae 时 刻 发 生 破 碎 . 在 71% 设 计 转 速 条 件 下 RI 起 始 来 源 
力 面 形成 的 ， 如 此 通道 3- 作 申 的 流 场 随时 间 做 周期 ” 出 现 了 泄漏 涡 破碎 的 参与 . 即 随 着 转速 的 提高 ，RI 
性 变化 . 起 始 的 过 程 发 生 馈 灾 , 进而 RI 的 参数 也 随 之 改变 . 
上 述 固定 通道 中 泄漏 流 的 周期 性 波动 正 与 六 
Zhang 等 0 和 Du 等 03 均 所 提 的 泄漏 流 自 诱 
导 非 定常 性 是 是 一 致 的 。 也 就 是 说 ， 在 本 文 的 研究 。 
对 象 中 53% 设 计 转 速 下 RI 的 产生 是 源 于 通道 内 叶  \ 
尖 泄 漏 流 自我 诱导 非 定常 性 的 周 向 传播 ， ©: 
3.4.2 71% 设 计 转 速 近 失速 工 况 流 场 分 析 AGS 
图 11 给 出 了 71% 设 计 转 速 数值 模拟 失速 边界 工 
况 五 通道 一 个 扰动 周期 泄漏 流 及 各 叶片 压力 面 静 压 
分 布 图 , 7 的 选取 为 压力 的 近 波 谷 到 近 波 谷 ， 如 图 
8(b) 所 示 。 图 中 红色 圆圈 依旧 示 出 叶片 3 压力 面 低 
压 区 的 周期 性 运动 。 除 了 与 53% 设 计 转 速 相同 的 尾 
缘 泄 漏 流 的 周期 性 运动 之 外 ,在 71% 设 计 转 速 压气 。 图 12 71% 设 计 转 速 数 值 模拟 失速 边界 工 况 t= (4/5)T 时 刻 通 


机 流 场 中 还 出 现 瞬时 的 主 泄漏 涡 破 碎 , 如 图 12 所 示 道 3-4 泄漏 流 线 图 
是 图 11(e) 中 通道 3-4 的 详细 泄漏 流 线 图 . Fig. 12 The TLF streamlines in passage 3-4 at t = (4/5)T for 


NS at 71% design rotor speed 


结合 图 11(e) 和 图 12 可 知 在 通道 3-4 中 叶片 5 
叶 顶 泄漏 流 发 生 第 二 次 泄漏 和 叶片 4 叶 顶 泄漏 流 的 
共同 作用 使 得 叶片 4 的 主 汇 漏 发 生 破碎 。 改变 转速 ”4 结 论 
其 实质 是 改变 了 叶片 的 载荷 ,叶片 叶 顶 载荷 的 增 大 
使 得 叶 顶 泄漏 流 强 度 增强 ,最终 使 得 主 泄漏 在 某 些 以 某 亚 音 轴 流 压气 机 转子 为 研究 对 象 ， 两 个 不 
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同 转速 的 试验 以 及 数值 模拟 研究 结果 结论 如 下 : 
1) 试验 测量 信号 分 析 表 明 ， 两 个 转速 条 件 均 
在 靠近 失速 边界 的 极 小 流量 范围 的 工 况 中 检测 到 非 
定常 性 RI。 在 近 失 速 稳定 工 况 ， 机 匣 壁 压力 信和 号 
频谱 存在 驼峰 特征 频率 ， 其 驼峰 峰值 处 频率 随 转速 
的 提高 相对 减 小 。 其 中 53% 设 计 转 速 驼峰 峰值 处 
频率 为 0.47BPF，71% 设 计 转 速 驼 峰 峰 值 处 频率 为 
0.34BPF. 
2) 数值 模拟 数值 监测 点 的 静 压 信号 分 析 表 明 ， 
两 个 转速 均 在 近 失 速 边 界 的 极 小 流量 范围 的 工 况 中 
叶 尖 区 存在 着 非 定常 性 。 其 在 相对 坐标 系 下 的 频率 
分 别 为 0.33BPF 和 0.26BPF, 在 绝对 坐标 系 下 的 频 
率 分 别 0.47BPF 和 0.34BPF, 这 与 试验 具有 一 致 性 ， 
说 明 五 通道 的 数值 模拟 成 功 的 预测 出 了 RI. 
3) 数值 模拟 流 场 分 析 说 明了 ， 叶 片 压 力 面 的 低 
压 区 周期 性 运动 影响 了 叶 顶 的 载荷 分 布 ， 从 而 使 得 
尾 缘 泄 漏 流 作 周期 性 的 变化 ， 尾 缘 泄 漏 流 冲击 相 邻 
四 厅 怀 加 轴 和 0 
终 形成 了 叶 尖 区 的 的 周 向 行 波 。 
4) 两 种 转速 下 , RI 均 源 于 泄漏 流 的 
定常 性 的 轴 向 传播 . 53% 设 计 转 速 情况 不 已 


71 儿 设计 转速 下 , 流 场 叶 尖 泄漏 流 自 FE 定常 性 
周期 性 运动 的 过 程 中 , 伴随 着 瞬 漏 涡 破碎 。 
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